import pandas as pd
import holoviews as hv

hv.extension('bokeh')

df_analytic = pd.read_csv('./Analytic.csv')
data = hv.Dataset(df_analytic)

analytic_phase_diagram = data.to(hv.HeatMap,['J','Delta'],'absmin_cos',groupby=['h'])
analytic_phase_diagram.opts(hv.opts.HeatMap(tools=['hover'],colorbar=True, width=600, height=600,toolbar='above',title='Analytic QPD',xrotation=90))
analytic_phase_diagram

h:

df_quspin = pd.read_csv('./experiments/quspin_experiment/quspin6.csv')

data2 = hv.Dataset(df_quspin)

numerical_phase_diagram = data2.to(hv.HeatMap,['J','Delta'],'M^2', groupby=['h'])
numerical_phase_diagram.opts(hv.opts.HeatMap(tools=['hover'],colorbar=True, width=600, height=600,toolbar='above',title='Numerical QPD N=6',xrotation=90))
numerical_phase_diagram

h:

df_bdg = pd.read_csv('./experiments/bdg_experiment/bdg_eigenvals_eigengraph6', converters={'index': str})

model = hv.Dataset(df)

bars = model.to(hv.Bars,'index',['ground','ground_J=0'],groupby=['J','Delta','N'])
bars.opts(hv.opts.Bars(tools=['hover'],width=700,xrotation=90,color='ground_J=0',cmap='Set3'))

J:

Delta:

N:

df = pd.read_csv('experiments/bdg_experiment/bdg_eigenvals_full_minvals6.csv', converters={'index': str})
data = hv.Dataset(df)
phase_diagram = data.to(hv.HeatMap,['J','Delta'],'M^2',groupby=[])
phase_diagram.opts(hv.opts.HeatMap(tools=['hover'],colorbar=True, width=600, height=600,toolbar='above',title='BDG QPD N=6',xrotation=90))
phase_diagram

df = pd.read_csv('experiments/bdg_experiment/bdg_eigenvals_full_minvals12.csv', converters={'index': str})

datapiv = df.pivot(index="Delta", columns="J", values="M^2")

phase_diagram = hv.Image((datapiv.columns,datapiv.index,datapiv.values))
phase_diagram.opts(hv.opts.Image(tools=['hover'],colorbar=True, width=600, height=600,toolbar='above',cmap='spectral',title='BDG QPD N=12',xlabel='J',ylabel='Delta',clabel='M^2'))
phase_diagram* hv.Text(0.5, 1, 'Ferromagnetism').opts(color='w')*hv.Arrow(0, 1.5, 'Antiferromagnetism', '<').opts(hv.opts.Arrow(text_color='w'))